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ABSTRACT 



Aims. To develop a method to determine distances to molecular clouds using JHK near-infrared photometry. 

Methods. The method is based on a technique that aids spectral classification of stars lying towards the fields containing the clouds 
into main sequence and giants. In this technique the observed (J - H) and {H - K s ) colours are dereddened simultaneously using trial 
values of A v and a normal interstellar extinction law. The best fit of the dereddened colours to the intrinsic colours giving a minimum 
value of^- 2 then yields the corresponding spectral type and A v for the star. The main sequence stars, thus classified, are then utilized 
in an A v versus distance plot to bracket the cloud distances. 

Results. We applied the method to four clouds, L1517, Chamaeleon I, Lupus 3 and NGC 7023 and estimated their distances as 
167 + 30, 151 ± 28, 157 ± 29 and 408 ± 76 pc respectively, which are in good agreement with the previous distance estimations 
available in the literature. 

Key words. ISM: dust extinction - ISM: clouds - Stars: distances - ISM: individual objects: L1517 -ISM: Chamaeleon I -ISM: 
individual objects: Lupus 3 - ISM: individual objects: NGC 7023 - infrared: ISM 
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The knowledge of distance to interstellar clouds which are ei- 
ther isolated or associated with large complexes, star forming 
or starless, is vitally important to provide observational con- 
straints to the properties of clouds and pre-main sequence (PMS) 
stars (e.g., Yun & Clemens 1990; Clemens, Yun & Heyer 1991; 
Kauffmann et al. 2008). However, distances to many of these 
clouds, especially those which are isolated, are highly uncertain 
(in some cases, by a factor of 2, see Hilton & Lahulla, 1995). 

In recent years, efforts have been made successfully to es- 
timate distances to large complexes using a number of meth- 
ods. For example, combining extinction maps from the Two 
Micron All Sky Survey (2MASS) with Hipparcos and Tycho par- 
allaxes, Lombardi et al. (2008) estimated accurate distances to 
the Ophiuchus and Lupus dark cloud complexes; the Hipparcos 
parallax measurements and B-band polarimetry were utilized 
by Alves & Franco (2007) to estimate the distance of the Pipe 
nebula. Recently, using multi-epoch very long baseline array 
(VLBA) observations of a number of young stellar objects 
(YSOs), Loinard et al. (2007) measured the trigonometric paral- 
laxes and hence distances to them with a precision of about 1-4% 
(Loinard et al. 2007, 2008). The distances to the parent clouds 
which harbour these YSOs are also constrained with similar pre- 
cision. 

However, it is difficult to establish distances to isolated 
clouds devoid of low-mass YSOs and primary indicators such 
as ionizing stars or reflection nebulae. Distances to these clouds 
can be assigned through their association (in position-velocity 
space) with nearby larger molecular clouds (e.g., Kauffmann 



et al. 2008) which have their distances determined through the 
above mentioned methods. The traditional ways of determining 
distances to isolated clouds utilized the star count method (Bok 
& Bok 1941) or Wolf diagrams (Wolf 1923) which is limited 
due to its dependence on the questionable extrapolation of lumi- 
nosity functions. Additional methods of assigning distances to 
small dark clouds entail bracketing the cloud distance by using 
spectroscopic distances to stars (Hobbs et al. 1986) and from the 
equivalent widths of interstellar Ca II H and K lines of stars close 
in front of and behind the cloud (Megier et al. 2005). However, 
spectroscopy of a large number of stars projected on to a cloud 
is extremely demanding in terms of telescope time. Because of 
the absence of nearby associated large molecular clouds or stel- 
lar spectroscopic data, distances of most of the isolated small 
clouds remain unknown. 

Distances to clouds can be determined by bracketing them 
using spectral classification of stars projected on to the cloud 
purely by photometry carried out in the Vilnius photometric sys- 
tem (Straizys, 1991; Straizys etal. 1992). f The distance and red- 
dening for individual stars can also be obtained using the very 
well calibrated Stromgren uvby — B intermediate-narrow- band 
photometric system. The physical parameters of stars thus ob- 
tained are utilized to bracket the cloud distances (e.g., Nielsen 
et al. 2000; Franco 2002). Another innovative method was de- 
veloped by Peterson & Clemens (1998) by identifying M dwarfs 
lying both in front of and behind the cloud. The reddening and 
spectral types of these M dwarfs were determined from photom- 
etry alone due to their conspicuous position in (V - T) versus 
(B-V) colour-colour (CC) diagram. The optical V, R, I and near- 
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IR (NIR) J,H &K 2MASS colours were utilized by Maheswar 
et al. (2004, 2006) to bracket cloud distances, again, by spec- 
tral classification of stars projected on to them from photometry 
alone. 

The enormous NIR data provided by the 2MASS 
(Kleinmann et al. 1994) were extensively used by several au- 
thors, mainly, (a) to discriminate between normal stars and stars 
with significant NIR emission from circumstellar material (e.g., 
Itoh et al. 1996) and (b) to map the extinction in dense molecu- 
lar clouds (Cambresy et al. 2002; Lombardi et al. 2006). Dutra 
et al. (2003) build an Ak extinction map of a 10° field towards 
the Galactic centre using 2MASS J and K s magnitudes. They 
utilized the upper giant branch of colour-magnitude diagrams 
and their dereddened mean locus built from previous studies of 
bulge fields to estimate the extinction. In this work we present 
a method by which distances to molecular clouds are estimated 
using (J-H), (H-K s ) colour indices of the main sequence stars 
in the spectral range from AO to K7. The technique is based 
on a methodology that provides spectral classification of stars 
projected on to the fields containing the cloud into main se- 
quence and giants by simultaneously dereddening the (J - H) 
and (H - K s ) observed colours using trial values of Ay and a 
normal interstellar extinction law. The main sequence stars are 
then utilized in an Ay vs. distance (d) plot to bracket the cloud 
distances. In ^2]we describe the method for finding the spectral 
type, Ay, and finally distances to the molecular clouds using NIR 
data from the 2MASS. The uncertainties involved in the deter- 
mination of the Ay values and the distances from the method are 
also discussed in this section. In §[3]we describe the criteria used 
to extract the data from the 2MASS database. In §|4] we apply 
our technique to four clouds, L1517, Chamaeleon I, Lupus and 
NGC 7023. We discuss the method in determining distances to 
dark molecular clouds in §|5] Finally, we conclude the paper by 
summarizing the results obtained in §|6] 

2. The Method 

The photometric distance d to a star is estimated using the dis- 
tance equation 

d(pc) = \Q( v - m v+5-Av)/5 _ \q(K-M k +5-A k )I5 ^) 

where V, My, Ay, K, Mk, and Ak are the apparent magni- 
tude, absolute magnitude and extinction in V and K filters re- 
spectively. Spectral type and luminosity class information are 
needed to determine both absolute magnitude of a star and the 
extinction suffered by it due to interstellar dust in the line of 
sight. Spectral type and luminosity class also determine vari- 
ous colours of a star. Therefore, in principle, one could deter- 
mine the spectral type and the luminosity of a star from its ob- 
served colours provided the extinction is zero. In practice, ob- 
served colours are reddened due to interstellar extinction which 
is wavelength-dependent. By assuming a value for Ay and an ex- 
tinction law one could estimate the colour excesses and correct 
the observed colours to find the intrinsic colours of the stars and 
hence their spectral types. 

2.1. The (J - H), (H - K s ) intrinsic colours 

In FigQ] we present the (J-H), (H-K s ) CC diagram. The solid 
line represents locations of unreddened main sequence stars. 
Usually the unreddened colours of different spectral types and 
luminosity classes in the 2MASS system are obtained from the 
Koornneef (1983) or Bessell & Brett (1988) using the transfor- 
mation equations given by Carpenter (2001). We produced the 



Table 1. The estimated (J - H)/„ and (H - K s )i„ intrinsic colours 
for main sequence stars. These colours were estimated using 
2MASS data of stars with known spectral types and distances 
from the Hipparcos parallaxes. Only those stars located within 
100 pc from the Sun were chosen. The N* gives the number of 
stars used to estimate the colours. 



Sp. 


M K 


(J - H) in 


(H — K s )i n 


(T (JH) n 


CT (HK)„ 


N* 


A0V 


0.65 


-0.026 


0.029 


0.030 


0.028 


50 


A1V 


1.00 


-0.014 


0.036 


0.037 


0.032 


45 


A2V 


1.17 


-0.004 


0.043 


0.044 


0.035 


43 


A3V 


1.35 


0.020 


0.048 


0.049 


0.032 


58 


A4V 


1.49 


0.042 


0.050 


0.044 


0.026 


23 


A5V 


1.60 


0.050 


0.052 


0.040 


0.029 


24 


A6V 


1.62 


0.058 


0.051 


0.037 


0.024 


12 


A7V 


1.63 


0.075 


0.054 


0.039 


0.031 


23 


A8V 


1.64 


0.084 


0.055 


0.047 


0.032 


8 


A9V 


1.72 


0.090 


0.059 


0.033 


0.026 


37 


F0V 


1.86 


0.119 


0.062 


0.042 


0.026 


105 


F1V 


2.00 


0.142 


0.062 


0.040 


0.026 


9 


F2V 


2.15 


0.162 


0.064 


0.035 


0.025 


141 


F3V 


2.26 


0.170 


0.065 


0.032 


0.029 


269 


F4V 


2.38 


0.175 


0.064 


0.037 


0.028 


28 


F5V 


2.49 


0.195 


0.068 


0.036 


0.027 


492 


F6V 


2.62 


0.211 


0.072 


0.036 


0.027 


390 


F7V 


2.75 


0.228 


0.069 


0.034 


0.032 


387 


F8V 


2.88 


0.235 


0.069 


0.038 


0.029 


308 


F9V 


3.04 


0.252 


0.080 


0.030 


0.028 


30 


GOV 


3.19 


0.257 


0.076 


0.033 


0.028 


396 


G1V 


3.24 


0.265 


0.078 


0.035 


0.028 


242 


G2V 


3.29 


0.276 


0.079 


0.036 


0.026 


320 


G3V 


3.34 


0.285 


0.083 


0.036 


0.029 


477 


G4V 


3.45 


0.302 


0.084 


0.032 


0.031 


8 


G5V 


3.56 


0.303 


0.085 


0.035 


0.030 


499 


G6V 


3.68 


0.325 


0.090 


0.040 


0.028 


188 


G8V 


3.90 


0.347 


0.093 


0.048 


0.030 


259 


G9V 


4.03 


0.363 


0.094 


0.074 


0.021 


9 


K0V 


4.15 


0.388 


0.098 


0.048 


0.031 


207 


K1V 


4.13 


0.406 


0.094 


0.046 


0.032 


93 


K2V 


4.15 


0.453 


0.115 


0.053 


0.028 


140 


K3V 


4.15 


0.473 


0.120 


0.046 


0.039 


137 


K4V 


4.22 


0.529 


0.130 


0.057 


0.037 


43 


K5V 


4.35 


0.555 


0.142 


0.051 


0.031 


108 


K7V 


4.83 


0.604 


0.160 


0.031 


0.035 


54 



intrinsic colours of the main sequence stars in the spectral range 
from AO to K7 in 2MASS system directly from the observa- 
tions by selecting stars from the All-sky Compiled Catalogue 
of 2.5 million stars (Kharchenko 2001) catalogue with known 
spectral types and which lie within 100 pc (estimated from their 
Hipparcos parallax measurements) from us. The mean values of 
(/ - H) and (H - K s ) colours produced are listed in columns 3 
and 4 of the Table Q~l respectively. The standard deviations and 
the number of stars used to evaluate the colours are given in 
columns 5, 6 and 7, respectively. The Mk values obtained from 
My and (V — K) from Cox (2000) are listed in column 2. The in- 
trinsic colors of the main sequence stars and giants in the spectral 
range K8 to M6 were converted from Bessell & Brett (1988) to 
the 2MASS system using the transformation equations given by 
Carpenter (2001). We interpolated the intrinsic colour indices of 
normal main sequence stars and giants to generate values with a 
0.5 spectral sub-type interval. 

The locations of the stars with various spectral types are 
identified and marked in FigJT] The dashed line is drawn paral- 
lel to the Rieke & Lebofsky (1985) interstellar reddening vector. 
An unreddened main sequence normal star would move paral- 
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Fig. 1. The (J - H) vs. (H - K s ) CC diagram drawn for stars 
towards a region chosen arbitrarily on the sky to illustrate the 
method described in §|2] The solid curve represents locations of 
unreddened main sequence stars. The reddening vector for an 
AOV type star drawn parallel to the Rieke & Lebofsky (1985) 
interstellar reddening vector is shown by the dashed line. The 
locations of the main sequence stars of different spectral types 
are marked with square symbols. The region to the right of the 
reddening vector is known as the NIR excess region and corre- 
sponds to the location of PMS sources. The dash-dot-dash line 
represents the loci of unreddened CTTSs (Meyer et al. 1997). 
The closed circles represent the observed colours and the arrows 
are drawn from the observed to the final colours obtained by the 
method for each star. 



lei to the reddening vector because of the extinction due to the 
interstellar dust along the line of sight. But as stars of spectral 
range AO to K7 move parallel to the reddening vector, it would 
be difficult to differentiate between reddened stars in the spectral 
range AO to K7 from unreddened M-type normal stars because 
the M-type stars lie across these vectors. However, for the stars 
lying below M-type stars loci, it is certain that they are reddened 
or unreddened normal stars in the spectral range AO to K7. The 
stars which occupy the region enclosed by the main sequence 
star loci (i.e., from A0-K7), reddening vector for an AO type star 
and below M-type loci were utilized in this work to determine 
distances to the dark clouds. 

2.2. Spectral types & A v from near-IR colours 

We used a procedure known as "photometric quantification", 
which is the procedure of inferring spectral type, and hence ab- 
solute magnitude and intrinsic colours of normal main sequence 
stars and giants from NIR photometry. In this method, stars with 
their (J -K s )< 0.750 were chosen from the fields containing the 
cloud. Then a set of dereddened colours [(J -H) and (H-K s ) ] 

1 This criterion allows us to avoid M-type stars from the analysis. 
Also, the classical T Tauri stars (CTTSs), often found to be associated 
the molecular clouds, occupy a well defined loci in NIR-CC diagram 
(Meyer et al. 1997) as shown in Fig. [TJ which intercepts with the main 



Fig. 2. The plot shows the relation between the uncertainties in 
(/ - H) and (H - K s ) colours of the sources in the four parts of 
the sky studied in this work. The Spearman rank-order correla- 
tion coefficient, r s , calculated for the regions are also given. The 
r s values indicate a strong correlation between the uncertainties 
of the two colours. The sources with the uncertainties < 0.035 
magnitude in J, H, & K s filters were considered. 

was produced for each star from their observed colours [(/ - H) 
and (H - K s )] by using trial values of Ay and the Rieke & 
Lebofsky (1985)0 extinction law in the equations 

(J-H) = (J-H)- 0.107 X Ay, (2) 

(H - K s ) = (H - K s ) - 0.063 X Ay, (3) 

the trial values of Ay were chosen in the range 0-10 magnitude 
with an interval of 0.01 magnitude though it is evident from 
Fig-fJJthat the maximum extinction that could be traced by this 

sequence loci at (J - K s )x 0.75. Thus the criterion of {J - K s )< 0.75 
would enable us to reject most of the CTTSs as well. 

2 The Rieke & Lebofsky (1985) law was derived towards the Galactic 
center in the Arizona-Johnson, not 2MASS, photometric system. But 
the extinction law derived towards the Galactic center in 2MASS sys- 
tem by Nishiyama et al. (2009) though consistent with the transformed 
Rieke & Lebofsky (1985) values, is different from those estimated by 
Indebetouw et al. (2005) derived towards massive star forming and 
"field" regions. It was suggested that the extinction law derived to- 
wards the Galactic center might be slightly different from those to- 
wards other sight lines. Since the maximum extinction that can be 
traced by our method is limited to ~ 4 magnitude, the regions stud- 
ied by us are in close resemblance to those by Indebetouw et al. (2005). 
However, the ratios, Aj/A K = 2.51 and A H /A K = 1.56 from the Rieke 
& Lebofsky (1985) values are in good agreement with the average val- 
ues of Aj/A K = 2.5 + 0.15 & A H /A K = 1.55 ± 0.08 obtained by 
Indebetouw et al. (2005) suggesting that the untransformed values of 
the Rieke & Lebofsky (1985) do represent the reddening law in the 
2MASS system. Similar suggestion was made by Alves et al. (1998) for 
their observations in CIT system. We used the values: Aj/A v = 0.282, 
A H /A V = 0.175, A K JA V = 0.112, following Cambresy et al. (2002) 
who found that the slope of the reddening vector measured in the 
2MASS colour-colour diagram is in better agreement with the Rieke 
& Lebofsky (1985) extinction law. 
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method is limited to Ay ~ 4 magnitude. The computed set of 
dereddened colour indices for a star were then compared with 
the intrinsic colour indices [(J - H)j„ & (H - K s )i n ] of nor- 
mal main sequence stars and giants, produced using the proce- 
dure discussed in the section |2~T1 The best fit of the dereddened 
colours to the intrinsic colours giving a minimum value of x 2 
then yielded the corresponding spectral type and Ay for the star. 
It is known that for a set of N data points with associated un- 
corrected uncertainties cr^., the best fit of predicted values yk.mod 
(from a model) to these points could be obtained by x 2 mini- 
mization defined as (e.g., Gould 2003) 



X 



z 

k=\ 



(yk - yk,mod) 



(4) 



But while the uncertainties in the magnitudes on /, H & K s fil- 
ters could be uncorrected, the uncertainties in the two colours 
(J - H) & (H - K s ) calculated from them with one filter (H) in 
common would be correlated. In Fig. [2] we show the correlation 
between the uncertainties in colours of stars from the parts of 
the sky (L1517, Chamaeleon I, Lupus 3 & NGC 7023) selected 
in this work. The Spearman rank-order correlation coefficient, 
r s , calculated for the sources in respective regions indicates a 
strong correlation between the uncertainties in the {J - H) and 
(H - K s ) colours. For correlated uncertainties, equation|4]is not 
valid and one has to use (e.g., Gould 2003) 

N N 

X 1 = XI Zj^* ~ ykjnod&kfal ~ yi,mod) (5) 
k=l 1=1 

where B = C _1 , the inverse of the colour covariance matrix. 
In the case of the three bands considered here, the colour co- 
variance matrix is given as (Lombardi & Alves 2001; Mafz- 
Apellaniz 2004) 



C 



J-H.H-K 



-<r H a- H + o~ K 



(6) 



with the off-diagonal elements having non-zero values unlike in 
the case of uncorrected uncertainties. The expected value of 
X 2 after it has been minimized is given by the difference be- 
tween the number of colours used and the number of parame- 
ters to be estimated (e.g., Gould 2003). Because we used two 
colours, (J - H) and (H - K s ), to estimate two parameters, spec- 
tral type and extinction, the expected values ofx 2 should be zero. 
However, since we used discrete values of Ay (0.01 magnitude) 
to deredden stars, the dereddened colours of the stars might not 
match exactly with the intrinsic colours of main sequence and 
giants to give x 2 values equal to zero but values <K 1 . We con- 
sidered only those solutions which gave a minimized^ 2 < 0.1 
for our analysis. 

The whole procedure is illustrated in the Fig[T]where we plot 
the NIR-CC diagram for stars chosen from an arbitrary loca- 
tion on the sky. The arrows are drawn from the observed data 
points to the corresponding dereddened colours estimated using 
our method. The maximum extinction values that can be mea- 
sured using the method are those for A0V type stars (« 4 mag- 
nitude). The extinction traced by stars falls as we move towards 
more late type stars. 

2.3. Distances to dark molecular clouds 

Once the spectral type and the Ay of the stars are known, we 
calculate their distances using the equation 

d (pc) = 10«.-*r+5-Ar)/5 (7) 



We excluded the stars classified as giants because of the large 
uncertainty in their absolute magnitudes. The distance at which 
the extinction of the stars showed a sudden increase in an Ay vs. 
d plot was considered as the distance to the cloud. 

2.4. Limitations and Uncertainties 

The uncertainty in the determination of distances to the clouds 
depends on the uncertainty in the estimation of distances to the 
individual stars which is given by 



+ <+<)xW/2-17) 2 



(pc) 



(8) 



where cr Ks is the photometric uncertainty in K s band, <x Mjf is the 
uncertainty in the estimation of the absolute magnitude or the 
spectral types and cta k is the uncertainty in the Ak (= 0.112 x 
Ay) estimated by the method. The major sources of uncertainty 
are, however, the photometric errors and the uncertainty due to 
a different reddening law which are statistical and systematic in 
nature, respectively. 

We allowed a maximum photometric uncertainty of 0.035 
magnitude for the sources in the J, H, & K s bands. The un- 
certainty in the absolute magnitudes or the spectral types de- 
termined from our method was evaluated by applying it to stars 
with known distances from the Hipparcos taken from the cata- 
logue by Kharchenko (2001). We used stars located within 100 
pc from the Sun so that the effect of the extinction could be ne- 
glected. Then, the equation [8] applied to the Hipparcos stars be- 
comes 



CT dH =^(cT 2 K +Cr 2 MK )x(dl2.\l) 2 



(9) 



In Fig. [3] we show histograms of the difference (d p - d/, pc) 
between the distances obtained using spectral types from our 
method, d p (pc) and those from the Hipparcos parallax mea- 
surements, dh (pc). FigureQa) shows the distribution when stars 
with <t„ h /jth < 0.2 (ratio of standard error to Hipparcos par- 
allax) from Hipparcos parallax measurements were selected. A 
Gaussian fit to the distribution shows a Icr deviation of ~ 18 
pc. The scatter in the distribution reduced to ~ 1 1 pc (Fig. Ob)) 
when stars with o~ 7rH / 71h < 0.05 from Hipparcos parallax mea- 
surements were selected. The majority (85%) of the latter stars 
were found to be located at distances < 60 pc. We noticed that 
the mean values of the distributions shown in both Fig.[3ja) and 
Fig-Etb) are shifted from the expected value of zero. Statistically 
the error in the mean, given by the standard deviation /y/N (N 
is the number of stars here), in both the distributions (a) and 
(b) should be less than 0.36 and 0.47 respectively. The negative 
shift in (a) could be due to the unaccounted binary stars which 
could result in the underestimation of distances using photom- 
etry. Adopting a Icr uncertainty of 18% of the distance for the 
Hipparcos stars in the Eq. [9] the uncertainty in the ct Mk is es- 
timated to be w 0.4 which corresponds to an uncertainty in the 
spectral type of about four subclasses Q. We assumed 0.4 as the 
uncertainty in the cr Mf , while calculating the uncertainty in the 
distances for all the stars. 

The Ay values obtained for the stars from our method are 
those which when used to deredden the observed colours (in 



3 The photometric errors considered in this work in the J, H, & K s 
magnitudes from the 2MASS database include the corrected band pho- 
tometric uncertainty, nightly photometric zero point uncertainty, and 
flat-fielding residual errors (Cutri et al. 2003). 

4 The change in M K with respect to the spectral types in the range 
A0V-K7V could be fitted with a function, M K a* 0.1 x S p. type + 0.98. 
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-100 -50 50 100 

d -d (pc) 

p h yr ' 

Fig. 3. (a) Histogram of the difference between the distances of 
the main sequence stars from our method and those from the 
Hipparcos parallax measurements with o-„ h Iji h < 0.2 are shown 
with a Gaussian fit to the data over-plotted, (b) same as in (a) but 
for main sequence stars with cr JlH /n H < 0.05. 



equations [2] and gave minimum values of x 1 - Therefore the 
Ay values obtained are equivalent to the average value of Ay 
calculated using the individual colours, i.e., Ay = [E(J - H) x 
9.35 + E(H - K s ) x 15.87]/2. Since the colours are found to be 
correlated, the expression for evaluating the uncertainty in Ay is 

cr(A v ) = J4J\o% + T-9 2 .cr 2 HKs + 2 x 37 X cov(JH, HK) (10) 

where, a 2 m = cr 2 , + <x^, cr^ 1K = °"h + °"i an ^ C0V (JH, HK) = 
r s x o-jho~hk s - Using the maximum photometric uncertainty of 
0.035 magnitude allowed in the individual bands and the highest 
value of 0.93 for r s (obtained towards L1517), the maximum 
uncertainty in the Ay was found to be ~ 0.6 magnitude. The 
uncertainty in the distances to the individual stars, equation|8] is 




10 100 1000 

distance (pc) 

Fig. 4. To illustrate the procedure used to determine distances to 
the clouds from an Ay vs. d plot as described in the 32. 51 



therefore dominated by the uncertainty in the determination of 
the spectral types or Mk values. 

The true spectral type and the Ay of the stars evaluated from 
our method directly depend on the reddening or extinction law 
that we adopt. Unlike in visual and UV wavelengths where the 
extinction law is known to vary along different sight lines in 
the Galaxy (Mathis 1990), the similarity in the wavelength de- 
pendence of extinction for very different line of sights passing 
through different environments like the giant H II region (RCW 
49) and the "field" region (Indebetouw et al. 2005) suggests that 
the extinction law is almost universal in the near-IR regime con- 
sidered here. Consequently, the observed colours can be cor- 
rected for the reddening to determine the spectral type of the 
stars using a standard extinction law. But any deviation of the 
reddening law from the standard extinction law towards the re- 
gion containing the clouds could introduce a systematic error in 
the distances estimated using our method. 

2.5. Determination of cloud distance from A v -distance plot 

The presence of an interstellar cloud should reveal itself by a 
sudden onset of extinction at a certain distance. If the cloud 
is isolated and is the first interstellar feature encountered, the 
excess should almost equal to zero in front of the cloud (eg., 
Straizys et al. 1992; Whittet et al. 1997; Knude & H0g 1998; 
Alves & Franco 2006; Lombardi et al. 2008). In this technique, 
the distance to the cloud is typically estimated from the first star 
that shows a significant reddening. The estimated distance is thus 
regarded as an upper limit due to the lack of information on the 
closeness of this background star to the cloud. In our method, be- 
cause we used relatively large number of stars projected towards 
the cloud, we often found an almost continuous distribution of 
stars with high extinction beyond a certain distance. 

In Fig. |4] we show a schematic figure of an Ay vs. d plot 
for an ideal situation along a line of sight through the Galactic 
plane. The extinction suffered by the stars increases as a func- 
tion of their distance (typically, ~ 1.8 mag kpc~ l is the value 
considered for the standard mean rate of diffuse-ISM extinction 
in the Galactic plane near the Sun, Whittet 1992) represented by 
filled circles. If we group the stars in distance bins and calculate 
the mean A y of the stars for each bin, they would follow the blue 



5 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



curve as shown in Fig. [4] The open circles represent the situation 
when a cloud of constant extinction is introduced at a certain dis- 
tance in the field-of-view. Now, the mean Ay of the stars in the 
bins behind the cloud would show a jump equal to the extinction 
introduced by the cloud (red curve). The distance at which the 
first sudden increase in the mean Ay occurs could be considered 
as the distance to the cloud (marked by the vertical dashed line 
in Fig.|4|. 

In real cases, we grouped the stars into distance bins of 
bin width = 0.09 x distance. The centers of each bin were sep- 
arated by the half of the bin width. Because there exist very few 
stars at smaller distances, the mean of the distances and the Ay 
of the stars in each bin were calculated by taking 1000 pc as the 
initial point and proceeding towards smaller distances (see Fig. 
IT7T ). The mean distance of the stars in the bin at which a sig- 
nificant drop in the mean of the extinction occurred was taken 
as the distance to the cloud and the average of the uncertainty 
in the distances of the stars in that bin was taken as the final 
uncertainty in our distance determination for the clouds. For all 
the clouds studied in this work, the vertical dashed lines in Ay 
vs. d plots, used to mark the cloud distances, were drawn at dis- 
tances deduced from the above procedure (see Fig.fTTIi. The error 
in the mean values of Ay were calculated using the expression, 
standard deviation /Va where N is the number of stars in each 
bin. 

3. The Data 

We extracted J, H and K s magnitudes of stars from the 2MASS 
All-Sky Catalog of Point Sources (Cutri et al. 2003) which sat- 
isfied the following criteria, 

1. Photometric uncertainty cr < 0.035 in all the three filters. 

2. Photometric quality flag of "AAA" in all the three filters, i.e., 
signal-to-noise ratio (SNR) > 10. 

The maximum allowed uncertainty of 0.035 magnitude in the 
photometry was chosen primarily as a trade-off between the min- 
imum scatter in the data points and the availability of enough 
stars in the background and the foreground of the clouds to 
bracket the cloud distances. One added advantage of choosing 
the data with high quality was that the limiting magnitude be- 
comes brighter and therefore, the selected stars would sample 
the lower density regions of the clouds where a normal redden- 
ing law is believed to be the most valid (Straizys et al. 1982; 
Kandori et al. 2003). 

4. Distances to the molecular clouds L1517, 
Chamaeleon I, Lupus 3 and NGC 7023 

We applied our method to four clouds, L1517, Chamaeleon I 
(Cha I), Lupus 3 (Lup 3) and NGC 7023, to determine dis- 
tances to them. We preferred these four regions as they are struc- 
turally less complex and have their distances estimated previ- 
ously. Also, these clouds harbour stars (inferred from associ- 
ated nebulosity) for which distances from the Hipparcos par- 
allax measurements are available (van den Ancker et al. 1998; 
Bertout et al. 1999) which would help constrain the parent cloud 
distances better. In Table [2] we list the distances to them com- 
piled from the literature. The central coordinates, total number of 
stars selected after applying all the selection criteria (cr < 0.035, 
SNR > 10, & (J - K s )< 0.75) and the number of stars classi- 
fied as dwarfs and that were used to determine distances to the 
individual clouds are listed in columns 2-5 of the Table[3]respec- 
tively. 



Table 2. Distances to the clouds L1517, Cha I, Lup 3 and NGC 
7023 compiled from the literature. 



Cloud Distances quoted Associated Ref. 

Identification in the literature (pc) star 

L 1517 144+^ AB Aur I 

Chamaeleon 1 115-215,158 HD 97048 2,3 

150 ± 15, 175+g 4, i 

Lupus 3 130-170,140 HR 5999 5,6,7 

100, 208^, 140 ± 2 8,9,10 

NGC 7023 440, 429^ 9 , 288 ± 25 HD 200775 11,1 

288 ± 25 12 



1-Bertout et al. (1999);2-Schwartz (1991); 3-Franco (1991); 4-Whittet 
et al. (1997); 5-Murphy et al. (1986); 6-Krautter (1991); 7-Hughes et 

al. (1993); 8-Knude & H0g (1998); 9-van den Ancker et al. (1998); 

10-de Zeeuw et al. (1999); 11-Viotti (1969); 12-Straizys et al. 1992. 
' The Hipparcos distance to the AB Aur. 

Table 3. The central coordinates, the number of stars selected 
from each field after applying all the selection criteria (cr < 
0.035, SNR > 10, & (J - K s )< 0.75) and the stars classified 
as dwarfs using our method and used to determine distances to 
L1517, Cha I, Lup 3 and NGC 7023 are given. 



Field 


a(J2000) 


,5(72000) 


Total 


Stars classified 


Id. 


O 


O 


stars 


as dwarfs 


L1517 










Fl 


74.166973 


+30.557486 


395 


292 


F2 


73.670855 


+31.053637 


339 


260 


F3 


73.555237 


+30.138617 


508 


402 


F4 


73.080485 


+30.648516 


299 


230 


Total 






1541 


1184 


Chamaeleon I 








Fl 


169.399281 


-76.719727 


521 


360 


F2 


165.432583 


-76.312943 


494 


352 


F3 


167.607235 


-77.711487 


235 


161 


F4 


163.267247 


-77.305313 


360 


260 


F5 


165.209741 


-78.675217 


502 


351 


F6 


160.921746 


-78.242210 


550 


389 


Total 






2662 


1873 


Lupus 3 










Fl 


242.405769 


-39.179062 


490 


335 


F2 


241.008065 


-39.471130 


814 


561 


F3 


243.129521 


-38.354233 


487 


308 


Total 






1791 


1204 


NGC 7023 








Fl 


315.431610 


+68.160049 


310 


230 


F2 


314.744713 


+67.160805 


394 


281 



4.1. L 1517 

In Fig.[5]we present the 2° x 2° high resolution extinction map0 
produced by Dobashi et al. (2005) of the region containing 
L1517. The contours are drawn at 0.5, 1.0, 1.5 and 2 magni- 
tude levels. The Herbig Ae star, AB Aur (Herbig 1960; The et al 
1994), associated with the reflection nebulosity, GN 04.52.5.02 
(Magakian 2003), is identified in Fig. [5] 

Some of the stars classified as main sequence by our method 
could actually be giants. Such erroneous classifications could 
lead to an underestimation of their distances which could cre- 
ate confusion with the real increase in the extinction due to the 



5 The extinction map, covering the entire region in the galactic lati- 
tude range \b\ < 40° derived using the optical database "Digitized Sky 
Survey I" and the traditional star-count technique, was produced in two 
angular resolutions of 6' and 18' (Dobashi et al. 2005). 
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Fig. 5. The 2° x 2° extinction map produced by Dobashi et al 
(2005) containing L1517 is shown with the fields F1-F4, each 
covering 40' x 40' area, marked and labeled. The contours are 
drawn at 0.5, 1.0, 1.5 and 2.0 magnitude levels. The stars used 
for estimating distance to L15 17 are represented by filled circles. 
The position of the Herbig Ae star, AB Aur, is also identified 
and labeled. The star with high extinction (Ay > 1 magnitude) 
at a distance closer than the estimated distance of the cloud is 
identified using the star symbols. 
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Fig. 6. The NIR-CC diagrams for the stars selected from the four 
fields containing L1517 are shown. The dots represent all the 
stars in a given field with photometric errors < 0.035 magnitude 
in J, H, & K s bands and open circles represent stars with their 
(J - K s ) < 0.75. 
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Fig. 7. The A y vs. d plot for all the stars obtained from the fields 
F1-F4 combined together towards L1517. The dashed vertical 
line is drawn at 167 pc inferred from the procedure described in 
the {32.51 The solid curve represents the increase in the extinc- 
tion towards the Galactic latitude of b — -7.98° as a function 
of distance produced from the expressions given by Bahcall & 
Soneira (1980). The error bars are not shown on all the stars for 
better clarity. 



presence of a cloud. This problem could be circumvented if we 
divide a large field, containing the cloud, into small sub-fields. 
While the rise in the extinction due to the presence of a cloud 
should occur almost at the same distance in all the fields, if the 
whole cloud is located at the same distance, the wrongly clas- 
sified stars in the sub-fields would show high extinction not at 
same but random distances. We divided the region containing 
L1517 into four fields as shown in Fig. [5] Each field is of 40'x40' 
in area. The central coordinates, the number of stars selected af- 
ter applying all the selection criteria (<r < 0.035, SNR > 10, 
& (J - K s )< 0.75) and the number of stars classified as dwarfs 
by our method from each field are listed in Table [3] In Fig. [6] 
we show the NIR-CC diagrams for the stars selected from Fl- 
F4. The stars with their photometric errors < 0.035 magnitude 
in /, H, & K s bands are shown as dots and among them the stars 
with (J - K s ) < 0.75 are identified by open circles. The stars 
classified as dwarfs using our method are identified by filled cir- 
cles (in white) in Fig. [5] It can be noticed that very few sources 
are selected towards the high extinction regions due to our data 
selection criteria discussed in section [3] 

In Fig. [7] we show the Ay vs. d plot for the stars from all 
the fields combined together. The solid curve shows the increase 
in the extinction towards the Galactic latitude b = -7.98° as 
a function of distance produced using the expressions given by 
Bahcall & Soneira (1980) (BS80, hereafter). Clearly, the stars 
with high extinction located beyond 167 pc, marked by the ver- 
tical dashed line (drawn using the procedure described in the 
{32.51 see Fig.fTTb. are distributed almost continuously. There ex- 
ist one star with Ay > 1 magnitude at distance closer than 167 
pc. The Ay vs. d plots for the stars from the individual fields 
F1-F4 are shown in the Fig. [8] The stars showing high extinction 
at 167 pc are located in the fields, Fl and F4, which contain the 
major parts of the cloud. The high extinction star closer than 167 
pc in the field F4 is identified using filled star symbol in the Fig. 
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Fig. 8. The Ay vs. d plots for the stars from the fields F1-F4 
towards L1517. The dashed vertical line is drawn at 167 pc in- 
ferred from the procedure described in the 32.51 The solid curve 
represents the increase in the extinction towards the Galactic lat- 
itude of b — -7.98° as a function of distance produced from the 
expressions given by Bahcall & Soneira (1980). 



[5] Because the star is projected within the cloud boundary, it is 
likely to be a misclassified source located behind the cloud. 

The distances to the molecular clouds in the Taurus-Auriga 
complex is considered as 140 pc based on the star count, photo- 
metric distances of reflecting nebulae and reddening versus dis- 
tance diagrams of the field stars towards the entire complex. But 
the most reliable estimation of the distance to L 15 17 is based on 
the Hipparcos parallax measurements of the two stars, AB Aur 
and SU Aur, found associated with the cloud. The Hipparcos 
distances to these stars are 144^ and 152^ respectively. The 
distance of 167 + 30 pc to L1517 determined from our method 
are found to be in good agreement with the distances of AB Aur 
and SU Aur within the uncertainties. 
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Fig. 9. The 5° x 5° extinction map produced by Dobashi et al 
(2005) containing Cha I is shown with the contours drawn at 
0.5, 1.0, 1.5 and 2.0 magnitude levels. The fields F1-F6 selected 
towards Cha I are marked and labeled in the figure. The stars 
classified as dwarfs using our method are identified with filled 
circles. Each field is of 1° x 1° in area. The locations of the stars 
HD 97300 & HD 97048 are identified and labeled. The stars 
with high extinction (A y > 1 magnitude) at distances closer than 
the estimated distance of the cloud is identified using the star 
symbols. 



: Cha 




distance (pc) 

Fig. 10. The Ay vs. d plot for the stars from all the fields (F1-F6) 
of Cha I combined together. The vertical dashed line is drawn 
at 151 pc inferred from the procedure described in the 32.51 The 
vertical dotted line is drawn to indicate the distance of a possible 
foreground dust layer at ~ 128 pc. The solid curve represents the 
increase in the extinction towards the Galactic latitude of b — 
-15.9° as a function of distance produced from the expressions 
given by BS80. For better clarity, the error bars are not shown 
on all the stars. 
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Fig. 11. The Ay vs. d plot for stars in the fields F1-F6 selected 
towards Cha I are shown here. The dashed vertical line is drawn 
151 pc inferred from the procedure described in the 32.51 The 
dotted vertical line is drawn a 128 pc. The solid curve represents 
the increase in the extinction towards the Galactic latitude of b — 
-15.9° as a function of distance produced from the expressions 
given by Bahcall & Soneira (1980). 



4.2. Chamaeleon 1 



The region containing the cloud Cha I is shown in the 5° x 5° 
extinction map produced by Dobashi et al (2005) in Fig. [9] The 
contours are drawn at 0.5, 1.0, 1.5 and 2.0 magnitude levels. 
Since the cloud is extended over a region of about 3° x 2° area, 
we divided the whole region into six sub-fields of 1° x 1° each 
as identified in Fig. [9] The central coordinates, total number of 
stars selected after applying all the selection criteria (<x < 0.035, 
SNR > 10, & (J - K s )< 0.75) and the number of stars classi- 
fied as dwarfs using our method in each field are listed in Table 
[3] Stars classified as dwarfs are identified in Fig. [9] using filled 
circles. Here too we find very few sources towards the high ex- 
tinction regions. 

In Fig. [10] we show the Ay vs. d plot for stars from all the 
fields combined together. The solid curve represents the increase 
in the extinction towards the Galactic latitude of b = -15.9° as 
a function of distance produced using the expressions given by 
BS80. Using the procedure described in the 32.51 we found a 
drop in the extinction at 151 pc (marked by the vertical dashed 
line in the Fig. [TOl as shown in the Fig. [T7] But the drop in the 
extinction at 151 pc was not as conspicuous as was observed to- 
wards L1517. Three stars with Ay > 1 were found to be located 
closer than 151 pc (identified in Fig. fusing filled star symbols). 
We noticed another drop in the extinction at 128 pc (see Fig. 
Wl\ . In Fig.QT] we show Ay vs. d plot for stars in the individual 
fields. The dotted line is drawn at 128 pc. The stars with ex- 
tinction significantly above the values expected from the BS80 
model seem to be present in all the fields at or beyond 15 1 pc un- 
like the distribution of the stars closer than 151 pc implying that 
the foreground dust layer at 128 pc is very patchy in distribution. 

Earlier attempts to determine distances to Cha I were made 
by Whittet et al. (1987), Franco (1991), Schwartz (1991) and 
Whittet et al. (1997). Schwartz (1991) reported a value of 115- 
215 pc, Franco (1991) assigned 158 pc while Whittet et al. 
(1997) deduced a most probable distance of 150 ± 15 pc based 
on the reddening vs. distance plot for stars towards Cha I. 
Independent distance estimates based on parallax measurements 
made by the Hipparcos satellite are available for two stars, HD 
97048 and HD 97300 (van den Ancker et al. 1998; Bertout et al. 
1999), which illuminate prominent reflection nebulae, Ced 111 
and Ced 112 respectively. The parallax distances of HD 97048 
and HD 97300 were estimated to be 175+^ and 188+^J respec- 
tively (Bertout et al. 1999). But the spectrophotometric distances 
to these stars were estimated to be 180 ± 20 and 152 ± 18 respec- 
tively (Whittet et al. 1997). These stars are located in F3 and 
Fl as identified in Fig. [9] Based on uvby/3 photometry of 1017 
stars, Corradi et al. (1997) showed the presence of a dust layer 
which is a part of a large scale distribution extended over a large 
area covering Chamaeleon, Musca and Coalsack at a distance 
of ~ 150 pc. In a subsequent work (Corradi et al. 2004) they 
showed the presence of gas as close as 60 pc from the Sun using 
interstellar Na I D absorption lines towards 63 B-type stars. The 
dust layer found in Fig. [10] and QT| at 128 pc could be a part of 
the forground dust layer shown by Corradi et al. (1997, 2004). 
The onset of extinction in Ay vs. d plots for both Fl and F3 are 
found to occur consistently at ~ 151 pc (Fig. ITTb. Based on our 
method, the present determination of 151 + 28 pc to Cha 1 is 
found to be in good agreement with the most probable distance 
estimated for the cloud by Whittet et al. (1997) and Corradi et 
al. (1997, 2004). 
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Fig. 12. The 4° x 4° extinction map produced by Dobashi et al 
(2005) containing Lup 3 is shown with the contours drawn at 
0.5, 1.0 , 1.5, 2.0, 2.5 and 3.0 magnitudes. The three fields F1-F3 
each of 1° x 1° area containing the cloud and the Herbig Ae/Be 
star, HR 5999, are identified and labeled in the figure. The stars 
with high extinction (Ay > 1 magnitude) at distances closer than 
the estimated distance of the cloud is identified using the star 
symbols. 



: Lup 3 




distance (pc) 

Fig. 13. The Ay vs. d plot for the stars from all the fields F1-F3 
of Lup 3 combined together. The dashed line is drawn at 157 pc 
inferred from the procedure described in the §12.51 The error bars 
are not shown on all the stars for better clarity. 



4.3. Lupus 3 

The region containing the Lup 3 cloud is shown in the 4° x 4° 
extinction map produced by Dobashi et al (2005) in Fig. [12] The 
contours are drawn at 0.5 to 3.0 with an interval of 0.5 mag- 
nitude. The three fields containing the cloud are identified in 
square boxes. The field, Fl contains Herbig Ae star, HR 5999 
(The et al 1994). The central coordinates, total number of stars 
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Fig. 14. The A y vs. d plots for stars from the individual fields Fl - 
F3 of Lup 3. The dashed line is drawn at 157 pc inferred from the 
procedure described in the §12.51 The dot-dashed line shows the 
increase in extinction towards the Galactic latitude of b — +9.4° 
as a function of distance estimated using expressions given by 
BS80. The source marked with an open triangle is identified with 
HD 143261, a giant (K1/K2III). 

selected after applying all the selection criteria (cr < 0.035, 
SNR> 10, & (J - K s )< 0.75) and the number of stars classified 
as dwarfs (identified in Fig. [12] with filled circles) are listed in 
Table[3] As seen in the cases of L1517 and Cha I, towards Lupus 
3 too, we found very few stars towards the dense parts of the 
cloud. 

In Fig. [13] we show the Ay vs. d plot for the stars from all the 
three fields combined. The solid curve represents the increase in 
the extinction towards the Galactic latitude of b — +9.4° as a 
function of distance produced using the expressions given by 
BS80. We found a sharp drop in the extinction at 157 pc, as 
shown in Fig. [17] estimated using the procedure described in the 
§12. 5 1 (marked with the vertical dashed line in the Fig. [131 . There 
are three stars with Ay > 1 at distances closer than 157 pc in Fig. 
[T3l In Fig. [14] we show the A y vs. d plot for the individual fields 
towards the Lup 3. The source located at 49 pc and Ay — 0.5 in 
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Fig. 15. The 3° x3° extinction map of the region containing NGC 
7023 is shown with the fields Fl and F2 ( of 1° x 1° each) iden- 
tified. The contours are drawn at 0.6, 1.0, 1.5, 2.0, 2.5 and 3.0 
magnitude. The arrow shows the location of the Herbig Ae/Be 
star HD 200775. The stars (only from Fl) with high extinction 
(Ay > 1 magnitude) at distances closer than the estimated dis- 
tance of the cloud is identified using the star symbols. 
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F2 is identified with HD 143261, a giant, classified as K1/K2III 
(Simbad database). Of the two stars with Ay > 1 in F3, the one 
located at 110 pc and Ay = 1.9 is identified with HD 145355 
classified as A1III (Simbad database). 

Most of the earlier studies have determined or assumed dis- 
tances in the range 130-170 pc either by considering Lupus to be 
associated with Scorpius-Centaurus association (Murphy et al. 
1986; Krautter 1991) or from spectroscopic parallaxes of early- 
type stars thought to be members of Scorpius-Centaurus asso- 
ciation (Hughes et al. 1993). But Knude & H0g (1998) based 
on extinction of stars located in the foreground and the back- 
ground of Lupus complex with their distances estimated from 
Hipparcos parallaxes, determined a closer distance of 100 pc to 
Lup 3. However, the Hipparcos parallax distance of HR 5999, 
illuminating the nebula GN 16.05.2 (Magakian 2003), is esti- 
mated to be 208^2 pc (van den Ancker et al. 1998) somewhat 
farther than most of the distance estimates. The best distance es- 
timation currently available for the Centaurus-Lupus subgroup 
of the Scorpius-Centaurus association is 140 + 2 pc derived by 
de Zeeuw et al. (1999) based on Hipparcos parallaxes. Using our 
method we determined a distance of 157 + 28 pc to Lup 3 cloud. 

4.4. NGC 7023 

The reflection nebula NGC 7023 is illuminated by the Herbig 
Be star HD 200775 (The et al. 1994). The 3° x 3° extinction map 
of the region containing NGC 7023 is shown in Fig. [15] The 
contours are drawn at 0.6, 1.0, 1.5, 2.0, 2.5 and 3.0 magnitude. 
Two fields of 1° x 1° were chosen towards the region. The central 
coordinates, total number of stars selected after applying all the 
selection criteria (cr < 0.035, SNR > 10, & {J - K s )< 0.75) and 



Fig. 16. The Ay vs. d plots for the stars in Fl and F2 towards 
NGC 7023. The dashed lines is drawn at 408 pc inferred from 
the the procedure discussed in the 32.51 The dotted lines indicate 
the locations of two additional dust layers possibly present at 
* 200 pc and =s 310 pc. The solid curve represents the increase 
in extinction towards the Galactic latitude of b = +14.1° as a 
function of distance estimated using expressions given by BS80. 
The error bars are not shown on all the stars for better clarity. 



the number of stars classified as dwarfs are listed in the Table [3] 
These stars are identified with filled circles in Fig. [TBI 

In Fig. [16] we present the Ay vs. d plot for stars from the 
fields Fl and F2. The solid curve shows the increase in the ex- 
tinction towards the Galactic latitude of b = +14.1° as a func- 
tion of distance produced using the expressions given by BS80. 
Clearly, a drop in the extinction is noticeable at two distances for 
Fl, at « 310 pc and at as 410 pc. In F2, though there exist a step 
like appearance again at distances » 310 and a; 410 pc, a number 
of stars show relatively high extinction at much closer distance 
of ss 200 pc, a distance where an increase in the Ay of > 1.5 
magnitude is seen in Fl also but for only two stars. From the Ay 
vs. d plot for stars in Fl and F2, we infer the presence of at least 
three layers of dust grains along the line of sight towards NGC 
7023. While the dust components at « 310 pc and * 410 pc are 
dominant towards Fl, the components at » 200 pc and * 310 pc 
are dominant towards F2. Using the procedure discussed in the 
§12.51 we found a significant drop in the mean value of Ay at 408 
pc (the distance marked by the vertical dashed line). We used the 
stars from Fl alone as the situation is much complex towards F2. 

The most widely quoted distances to LDN 1167/1174 or 
NGC 7023 and LDN 1147/1158 (clouds located ~ 2° west of 
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NGC 7023) groups are 288 ± 25 and 325 ± 13 pc respectively 
(Straizys et al. 1992). They used the Vilnius photometric system 
to classify the stars and also to get the interstellar reddening. 
They obtained a distance of 325 + 13 pc to LDN 1 147/1 158 by 
taking an average of 10 stars showing A y > 0.45 which were dis- 
tributed in the range 240-380 pc. For LDN 1 167/1 174 or NGC 
7023 group, they assigned a distance of 288 ±25 pc again by tak- 
ing the average distance of 4 considerably reddened stars. They 
preferred a distance of 275 pc to HD 200775 by assuming it to 
be a B3Ve star. 

The Hipparcos parallax distance to HD 200775 is estimated 
to be 429+^j 6 pc (van den Ancker at al. 1998; Bertout et al. 1999). 
The sharp rise in the extinction in Fig. [16] at 408 pc could be due 
to the association of the cloud with HD 200775 and the sharp rise 
in Ay at 305 pc could be due to a foreground dust component. 
Kun (1998) showed the presence of dust components at three 
characteristic distances: 200, 300, and 450 pc, based on cumula- 
tive distribution of field star distance moduli in Wolf diagrams. 
We found an additional component of dust at 170 pc towards the 
western parts of NGC 7023 especially towards LDN 1171 and 
LDN 1 147/1 158 group. The results will be presented in a forth- 
coming paper. Using a total of 230 stars, classified as dwarfs, we 
determined a distance of 408 ± 76 pc to the NGC 7023 cloud. 

5. Discussion 

The current determination of distances to L1517, Cha I, Lupus 
3 and NGC 7023 is based on the NIR photometry of relatively 
large number of stars, 1184, 1873, 1204 and 230, respectively 
(Table [3]), projected onto them. Distances estimated using our 
method to these clouds are found to be in good agreement with 
most of the previous estimations within the error (see Table[2j. 

Based on our method and the NIR data from the 2MASS, 
one could determine distances to any cloud located within ~ 500 
pc from the Sun (depending upon the Galactic latitude of the 
clouds). However, the complex nature of the interstellar medium 
is clearly evident in the Ay vs. d plots of Cha 1 and NGC 7023, 
as shown in Figures [10] & [16] respectively. There exist, clearly, 
at least two layers of dust components in the direction towards 
these two regions. If the first layer of dust component happens to 
be denser but uniform in distribution, it would reduce the proba- 
bility of finding the second dust layer. However, if the first layer 
is denser but patchy in distribution, we get sources with high ex- 
tinction at shorter distances before we could notice a wall of stars 
with high extinction due to the second layer of dust component. 
In that case, differentiating stars with spurious spectral classifi- 
cation from the actual rise in extinction due to the presence of 
dust components would become difficult. Dividing a larger field 
into sub-fields could help us decipher the Ay vs. d plots of a 
cloud. It is indeed noted that the presence of dominant dust com- 
ponents is revealed in Ay vs. d plots for different fields towards 
Cha 1 (Fig. [10| and NGC 7023 (Fig. [I6]l where the effects are 
most conspicuous. 

The size of the sub-fields required towards a cloud should be 
decided on the basis of its galactic location. We noticed that in a 
given field, our method consistently assigns as 70% of the total 
selected stars, after applying all the selection criteria (cr < 0.035, 
SNR > 10, & (J - K s )< 0.75), as dwarfs. Our experience from 
the application of the method to the four clouds in the previ- 
ous sections showed that we require at least 100-200 dwarfs to 
discern the presence of a sudden rise (or drop) in the extinc- 
tion in Ay vs. d plot of a field. The number of stars decreases 
with the Galactic latitude, thus, requiring a larger field as we go 
towards relatively high galactic latitudes. At higher galactic lat- 



< 



1 .A 
1 .2 
1 .0 

0.8 

CD 

E 0.6 

0.4 

0.2 
1 .2 

1 .0 

1 0.8 

g 0.6 
E 

0.4 

0.2 
1 .4 
1 .2 

> 1 .0 

c 

- 0.8 

\ 0.6 

0.4 

0.2 
1 .0 



0.8 - 



0.6 - 



CD 

E 0.4 




d= 1 5 1 pc 
Cha I 





d=157 pc 
Lup 3 




0.2 



100 200 300 400 500 600 
distance (pc) 

Fig. 17. The mean values of Ay vs. the mean values of distance 
plot for L1517, Cha I, Lupus 3 and NGC 7023 produced using 
the the procedure discussed in the 32.5l to determine distances to 
the clouds. The error bars on the mean A y values were calculated 
using the expression, standard deviation /y/N, where N is the 
number of stars in each bin. 



itudes, however, the confusion due to multiple dust components 
in the line of sight would be reduced as the majority of the clouds 
tend to occupy regions closer to the Galactic plane. 

For small and isolated clouds, it would be difficult to divide 
the field containing the cloud into sub-fields with sufficient num- 
ber of stars projected on to them and therefore the determination 
of distances made correspondingly hard. In such cases, clouds 
which are spatially closer and show similar radial velocities are 
to be chosen. Clouds located in the close proximity with similar 
radial velocities are believed to be at similar distances. 



6. Conclusions 

We present a method to determine distances to molecular clouds 
using the prodigious NIR data provided by the 2MASS. The 
method involves the following steps: 
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1. Extract /, H, K s magnitudes with photometric uncertainty 
< 0.035 and SNR > 10 of stars projected onto the fields 
containing the cloud from the 2MASS database. 

2. Select stars with their (/ - K s )< 0.75 to eliminate M-type 
stars from the analysis as unreddened M-type stars located 
across the reddening vectors of A0-K7 dwarfs make it dif- 
ficult to differentiate the reddened A0-K7 dwarfs from the 
unreddened M-type stars. 

3. A set of dereddened colours for each star is produced from 
their observed colours by using a range of trial values of Ay 
(0-10 mag) and the Rieke & Lebofsky (1985) extinction law 
in equations [2] and [3] 

4. Computed sets of dereddened colours of a star are then com- 
pared with the intrinsic colours of the normal main sequence 
stars. The intrinsic colours of the main sequence stars in the 
spectral range A0-K7 are produced from the 2MASS data of 
stars with known spectral types and those located within 100 
pc, calculated from the Hipparcos parallaxes. The best match 
giving a minimum value of y 1 as defined in the equation [5] 
then yields the spectral type and Ay corresponding to that 
intrinsic colour. 

5. Once the spectral types and A y values of the stars are known, 
their distances are estimated using the distance equation Q] 
Only those stars which are classified as dwarfs are consid- 
ered for the determination of distances as the absolute mag- 
nitudes of giants are highly uncertain. 

We applied the method first to main sequence stars which are 
located at distances < 100 pc (from Hipparcos parallaxes) and 
have their spectral types known. We showed that our method can 
determine the spectral types of the stars with a lcr uncertainty of 
about four subclasses which propagates into the uncertainty in 
the distances of the stars by < 20%. By applying the method 
to four clouds, L1517, Chamaeleon I, Lupus 3 and NGC 7023 
which harbour nebulous stars with the Hipparcos parallax dis- 
tances known, we determined a distance of 167 ± 30, 151 + 28, 
157 + 29 and 408 + 76 pc, respectively. These are found to be 
in good agreement with the most accurate distances available for 
them in the literature. We found additional components of dust 
layers at » 100 pc towards Chamaeleon I and at * 200 and ^ 310 
pc towards NGC 7023. In the future, we will extend the applica- 
tion of this method to many more molecular clouds to determine 
their distances using the 2MASS NIR photometry. 

Acknowledgments 

We thank the referee, Dr. Franco, G. A. P, for his useful com- 
ments on the work. This publication makes use of data products 
from the Two Micron All Sky Survey, which is a joint project 
of the University of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Technology, funded 
by the National Aeronautics and Space Administration and the 
National Science Foundation. This research has also made use of 
the SIMBAD database, operated at CDS, Strasbourg, France. S. 
D. is supported by the MAGNET project of the ANR (France). 

REFERENCES 

Alves, Joao., Lada, Charles J., Lada, Elizabeth A., Kenyon, 

Scott J., Phelps, Randy., 1998, ApJ, 506, 292 

Alves, F. O., Franco, G. A. P., 2007, A&A, 470, 597 

Bahcall J. N., Soneira R. M., 1980, ApJS, 44, 73 

Bertout, C, Robichon, N., Arenou, F, 1999,A&A, 352, 574 

Bessell, M. S., Brett, J. M., 1988, PASP, 100, 1134 



Bok B. J., Bok P. F, 1941, The Milky Way. Harvard Univ. Press, 
Cambridge , MA 

Cambresy, L., Beichman, C. A., Jarrett, T. H., Cutri, R. M., 

2002, AJ, 123, 2559 

Carpenter J. M., 2001, AJ, 121, 2851 

Clemens, Dan P., Yun, Joao Lin., Heyer, Mark H, 1991, ApJS, 
75, 877 

Corradi, W. J. B., Franco, G. A. P., Knude, J., 2004, MNRAS, 
347, 1065 

Corradi, W. J. B., Franco, G. A. P., Knude, J., 1997, A&A, 326, 
1215 

Cox, Arthur N., Allen's astrophysical quantities, 4th ed. 
Publisher: New York: AIP Press; Springer, 2000. Edited by 
Arthur N. Cox. ISBN: 0387987460 

Cutri R. M. et al., 2000, 2MASS All-Sky Catalog of Point 
Sources. NASA/IPAC Infrared Science Archive. 
Dobashi, Kazuhito., Uehara, Hayato., Kandori, Ryo., et al., 
2005, PASJ, 57S 

Dutra, C. M., Santiago, B. X., Bica, E. L. D., Barbuy, B., 2003, 
MNRAS, 338, 253 

de Zeeuw, P. T, Hoogerwerf, R., de Bruijne, J. H. J., Brown, A. 

G. A., Blaauw, A., 1999, AJ, 1 17, 354 

Franco, G. A. P., 1991, A&A, 251, 581 

Franco, G. A. P., 2002, MNRAS, 331, 474 

Gould, Andrew., 2003, astro.ph. 10577 

Herbig, George H, 1960, ApJS, 4, 337 

Hilton, J., Lahulla, J. F, 1995, A&AS, 113, 325 

Hobbs, L. M., Blitz, L., Magnani, L., 1986, ApJ, 306L, 109 

Hughes, Joanne, Hartigan, Patrick, Clampitt, Lori, 1993, AJ, 

105,571 

Indebetouw, R., Mathis, J. S., Babler, B. L., et al., 2005, ApJ, 
619, 931 

Itoh, Yoichi, Tamura, Motohide, Gatley, Ian., 1996, ApJ, 465L, 
129 

Kandori, Ryo., Dobashi, Kazuhito., Uehara, Hayato., Sato, 
Fumio., Yanagisawa, Kenshi., 2003, AJ, 126, 1888 
Kauffmann, J., Bertoldi, F, Bourke, T. L., Evans, N. J., II, Lee, 
C. W., 2008, A&A, 487, 993 

Kharchenko, N. V., Kinematika i Fizika Nebesnykh Tel, vol. 17, 
no. 5, p. 409-423 

Kleinmann, S. G, Lysaght, M. G, Pughe, W. L., et al. 1994, 
Ap&SS, 217, 11 

Knude, J., H0g, E., 1998, A&A, 338, 897 
Koornneef, J., 1983, A&A, 128, 84 

Krautter, J., Low Mass Star Formation in Southern Molecular 
Clouds, ESO Scientific Report. Edited by Bo Reipurth. 
Garching: European Southern Observatory (ESO), 1992., p. 127. 
Kun, M., 1998, ApJS, 115,59 

Loinard, Laurent, Torres, Rosa M., Mioduszewski, Amy J., 
Rodriguez, Luis F, Gonzalez-Lopezlira, Rosa A., Lachaume, 
Regis, Vazquez, Virgilio, Gonzalez, Erandy., 2007, ApJ, 671, 
546 

Loinard, Laurent, Torres, Rosa M., Mioduszewski, Amy J., 
Rodriguez, Luis F, 2008, ApJ, 675L, 29 
Lombardi, M., Alves, J., Lada, C. J., 2006, A&A, 454, 781 
Lombardi, M., Lada, C. J., Alves, J., 2008, A&A, 480, 785 
Magakian T. Y, 2003, A&A, 399, 141 

Maheswar G, Manoj P., Bhatt H. C, 2004, MNRAS, 355, 1272 
Maheswar, G, Bhatt, H. C, 2006, MNRAS, 369, 1822 
Mafz-Apellaniz, Jesus., 2004, PASP, 1 16, 859 
Mathis, John S., 1990, ARA&A, 28, 37 

Megier, A., Strobel, A., Bondar, A., Musaev, F. A., Han, Inwoo, 
KreAowski, J., Galazutdinov, G. A., 2005, ApJ, 634, 451 
Meyer, Michael R., Calvet, Nuria., Hillenbrand, Lynne A., 



13 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



1997, AJ, 114, 288 

Murphy, D. C, Cohen, R., May, J., 1986, A&A, 167, 234 
Nielsen, A. S., Jnch-Srensen, H., Knude, J., 2000, A&A, 358, 
1077 

Nishiyama, Shogo., Tamura, Motohide., Hatano, Hirofumi., et 
al., 2009, ApJ, 696, 1407 

Peterson, Dawn E., Clemens, Dan P., 1998, AJ, 1 16, 881 
Rieke, G. H., Lebofsky, M. J., 1985, ApJ, 288, 618 
Schwartz, R. D., 1991, in Low mass star formation in southern 
molecular clouds, ed. B. Reipurth, ESO Scientific Report no. 
11, p. 93 

Straizys, V., Wisniewski, W. Z., Lebofsky, M. J., 1982, Ap&SS, 
85,271 

Straizys, V., 1991, ppag., proc, 341 

Straizys, V., Cernis, K., Kazlauskas, A., Meistas, E., 1992 
BaltA, 1, 149 

The, P. S., de Winter, D., Perez, M. R., 1994, A&AS, 104, 315 
Torres, Rosa M., Loinard, Laurent., Mioduszewski, Amy J., 
Rodrguez, Luis E, 2007, ApJ, 671, 1813 

van den Ancker, M. E., de Winter, D., Tjin A Djie, H. R. E., 

1998, A&A, 330, 145 

Viotti, R., 1969, MmSAI, 40, 75 

Whittet, D. C. B., Kirrane, T. M., Kilkenny, D., Oates, A. P., 
Watson, F. G., King, D. J., 1987, MNRAS, 224, 497 
Whittet, D. C. B. 1992, Dust in the Galactic Environment 
(Bristol: IOP) 

Whittet, D. C. B., Prusti, T, Franco, G. A. P., Gerakines, P. A., 
Kilkenny, D., Larson, K. A., Wesselius, P. R., 1997, A&A, 327, 
1194 

Wolf M., 1923, Astron. Nachr., 219, 109 

Yun, Joao Lin., Clemens, Dan P., 1990, ApJ, 365L, 73 



14 



